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ABSTRACT 


Several models of ocean surface wave attenuation were employed in 
making hindcasts for comparison with spectra obtained from a wave 
recorder located at Point Sur, California. Quantitative decay models 
were reduced to common variables. The French Meteorological Service 
DSA II generation model was used. 

The attenuation model which best fits the data obtained is 
n = exp (-20 r t), where n is a ratio of energy at the start of 
decay to that remaining afterwards, T is wave period in seconds, 
and t is decay duration in hours. The period range considered was 


8.5 to 23.5 seconds. 
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I. INTRODUCTION 


Neither the magnitudes nor the mechanisms of ocean wave attenu- 
ation have been well determined for ocean surface waves with periods 
of 8.5 to 23.5 seconds. The standard wave forecasting method in the 
United States, for nearly 15 years, has been the Pierson, Neumann, 

and James [1955] method. The 'PNJ'’ technique allows for angular 
“spreading but fails to consider mechanisms of attenuation. 

Wave attenuation studies have been conducted using various 
methods of attack. Theoretical and experimental approaches lead to 
a variety of proposed mechanisms for the damping of wave energy. 
Wave forecasters have devised empirical methods of estimating loss 
of wave energy during propagation from a specified fetch area to 
an observation point. 

The French DSA series of forecasting methods attempted to 
improve on the limited fetch methods of 'PNJ’ by considering wave 
energy approaching the forecast point from several directions. The 
DSA method number two (DSA II) [Gelci, Cazale, and Vassal 1957] 
employs an attenuation model which allows for decay as a function 
of wave period, decay distance, and cross-wind condition. This 
method was developed in the North Atlantic Ocean. Field use in 
the Pacific Ocean provided an opportunity to evaluate DSA II under 
different geographical conditions. 

The objectives of this thesis were two-fold. First, to test DSA 
II in a region in the Northeast Pacific extending roughly 2000 miles 


from Point Sur, California and, second, to attempt to evaluate some 


de 


of the proposed attenuation models. The method employed was compari- 
son of hindcasts (made using DSA II) with analyzed wave records from 


Point Sur, California for January 1948. 
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II. ATTENUATION THEORIES AND MODELS 


When ocean waves propagate out of a region of active generation 
a period of attenuation commences. It is generally believed that 
longer waves decrease in height more slowly than their shorter 
counterparts. The following mechanisms and experiments provide 


-insight into the processes which may be involved. 


A. PROPOSED MECHANISMS OF DECAY 
l. Molecular Viscosity 
Viscous dissipation has been suggested. However, it has 
been noted theoretically that viscous amplitude decay time for ocean 
waves varies as the square of wave length [Kinsman 1965]. Except for 
very minute ocean waves, such as parasitic capillaries, viscous decay 
is several orders of magnitude too small to explain observed decay. 
2. Wave Breaking 
Breaking of waves is a method of converting wave energy into 
turbulence. [It is the most important mechanism in the generating 
area [Barnett 1966] and on the final beach but is rare elsewhere 
[Cartwright 1967]. Wave-wave interaction, as discussed below 
(II. A. 6.), may lead to breaking. 
3. External Turbulence 
Groen [1954] and Bowden [1950] ascribed wave damping primarily 
to externally caused turbulence. This turbulence is produced inde- 
pendently of the wave. An artificial "eddy vicsosity" may be used 
to deal with this mechanism as long as the size of turbulent eddies 


is small compared to the length of the waves [Barber and Ursell 1948]. 
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Bowden assumed monochromatic significant waves. Groen dealt with 

wave spectra to arrive at an energy decay expression of the following 
form: mn = exp - Geteoe Here n is the ratio of energy at the start of 
a region of decay to that present at the end of the region of length 

x Giaraie ib, mene deep water wave length in meters; a is Groen's 
attenuation coefficient in units of arse The energies were in 

a band centered at wave length L. 

Groen [1954, p. 7] stated that turbulence is frequency selective 
and attenuates the shorter period waves more rapidly. Thus it is 
better to employ wave spectra rather than assume, as a representative 
system, only the significant waves present. 

4. Air Resistance 

Darbyshire [1957] examined air resistance and concluded that 
it did not provide a dominant mechanism of decay since he found 
actual attenuation to be independent of wind speed for a following 
wind. 

Phillips [1966, p. 148] concluded that the condition of 9 
greater than 90° results ina very weak coupling between wind and 
swell which leads to an exponential form of decay. Here, 9 is the 
angle between the wind direction and the direction of energy propa- 
gation under consideration. He suggested that the seen ieee iS 
proportional to COS @Q. 

59. Parasitic Capillaries 

Phillips [1966, p. 134] discussed capillary waves which 
tend to form as standing waves in advance of the sharp crests of 


short gravity waves due to local surface tension effects. The 


capillary waves are small enough to cause viscous attenuation to 
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be significant. This mechanism is important in the case of steep 
gravity waves of length 5 to 30 cm. Therefore it is significant in 
small lakes but probably not in the ocean where most energy is con- 
centrated in lower frequencies. 

6. Wave-Wave Interaction 

Barnett [1966] suggested that, once the wind has ceased its 
active generation, the spectrum looses energy mainly through wave- 
wave interactions. Gelci and Cazale [1962, p. 30] have developed a 
model (DSA V) to consider losses due to the interaction between wave 
trains. DSA II does met consider this effect. Yet the decay 
observed in DSA V is quite similar to that resulting in DSA II. 

Gelci (personal correspondence) reports that wave-wave interaction 
does not significantly affect waves with a period of 14 seconds or 
reactor. 

The details of the wave-wave interaction mechanism are 
complex, but the effect is the transfer of energy from three active 
wave components to a fourth, passive component. The magnitude of 
the effect appears to be comparable to wave generation processes 
[Barnett 1966]. For a typical fully developed spectrum, energy is 
—— from a broad mid-frequency band to both higher and lower 
frequencies. The energy transferred to higher frequencies could be 
lost via breaking and viscosity. 

/. Scattering Due to Turbulence 

Phillips [1961] proposed that some ocean turbulence is 
generated by gravity waves. This is a type of turbulence not included 
in the externally caused turbulence (see II. A. 3.). Phillips [1959] 


developed a theory to explain attenuation and directional scattering 
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of surface gravity waves which appears to result from interaction 
between wave turbulence and wave motion. His conclusion was that 
scattering would predominate over viscous attenuation for wave 
lengths greater than about 3 meters (periods of 1.4 sec). It is not 
certain that this would be a significant effect since viscous effects 
themselves are small for the 8.5 to 23.5 second periods covered by 
this thesis. Phillips’ theory results in the following model for 


which values are tabulated in Appendix C: n = exp -(0.0087 pos oie 


B. EXPERIMENTS 

The first serious ‘application of modern analysis techniques to 
ocean waves occurred in the early 1940's [Cartwright 1967]. It was 
discovered that random fluctuations in sea level from wave Py 32): 
can be resolved with spectral analysis into the sum of a number of 
quasi-, »>ricodicoscillations of various frequencies or periods. Most 
of the wave measurements made during the next ten years were used to 
develop wave generation models and parameters for energy-density 
Specuera. 

1. San Clemente Study 

During the 1950's, Munk, Miller, Snodgrass, and Barber [1963] 

examined wave energy received at San Clemente, California from storm 
waves generated in distant regions of the North Pacific and in the 
Southern Ocean. They detected significant amounts of energy only 
at frequencies below 5.0 x 107° hz (periods larger than 20 sec). 
Wave energy at higher frequencies was heavily attenuated somewhere 
along the great circle path [Cartwright 1967]. This result was used 
in justification of hindcasting waves on the basis of surface winds 


within 2000 nm of Point Sur for this thesis. 


16 


2. New Zealand to Alaska Attenuation Study 


Munk and his colleagues set up one of the most ambitious 
series of measurements in the history of oceanography [Snodgrass 
et al. 1966]. During the summer of 1963 they tracked swell along a 
great circle path from New Zealand to Alaska using six observation 
stations. At frequencies above 7.0 x 107? hz (period 14.3 sec) there 
was a tendency for spectral energy to decrease at a rate of 1/10 db 
per degree along the great circle path. This is equivalent to the 
equation: 

n = exp [-(0.011 cia — c &\8 
Comparative values are presented in Appendix C. Below 7.2 x an hz 
the energy either remained constant or increased. Cartwright [1967] 
indicated that the major result of this experiment was that attenu- 
ation above 5.0 x ie hz (period 20 sec) in the first 20° or so from 
the generation area was much greater than that observed over the rest 
of the path to Alaska. 
3. Small Amplitude Deep Water Wave Study 

Carswell [1968], of the Naval Postgraduate School, looked at 
very small waves (frequencies 2 to 2.4 hz) in a lake. Following a 
suggestion of Professor J. B. Wickham, Carswell used his data to 
Tse er a exponential decay function. An N of three gave the best 

5 

fit and resulted in the following: n = exp | - (5 SO Fe x |. Here, 
f is frequency in hertz; X is decay distance in nautical miles. This 
expression was modified to the following equivalent form for evalu- 


ation in this thesis: 


S e -2 
n = exp | -(7.575 =) AL |. 


t/ 


Lee Spectro-Angular Attenuation 


DSA II has its own attenuation function, which was developed 
for use in the North Atlantic and is applicable where the angle 
between the wind direction in the decay region and the direction of 
energy propagation (@) exceeds 60° or wind velocity (W) is not greater 


than ten#knotsm Sbhis *tunct ion ws: 


n = exp [- L100 + 0.30816 (95) iia (Ge -45° psec” ea E 7 |, 
L 


When 9 does not exceed 60° and W is greater than ten knots DSA II 
assumes no decay. Gelci (personal correspondence) points out that 
this is nearly true for an isolated component and DSA II neglected 
"interferential decay" due to multiple wave trains. The case of @ 
no greater than 60° and W in excess of 10 knots is taken care of by 
either maintaining steady-state or generation as applicable for the 
wave period under consideration. 
5. Eddy Viscosity Application 
Groen and Dorrestein [1950] evaluated the attenuation 
coefficient (a), in the equation previously presented (II. A. 3.), 
1/6 


to get an empirical value (8.5 x 107" meter ). Their equation 


has® been modified tesaemore convenient f ormmhonmmthiis sehesis: 
ii ="exp | -(1.92 x en are eo) t pil) 


6. Laboratory Underwater Turbulence Study 


Paquin [1968] studied the effect of grid-generated underwater 
turbulence on small waves (1.2 to 12.3 hz) in a wave tank. His 


results showed that turbulent attenuation would have the following 


3a 
form: n = exp 2 RE] where K is a constant (].26 sec eee 
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X is a decay distance (in units compatible with g); g is gravitational 
acceleration; and f is frequency in hz. This was converted to vari- 


ables convenient to DSA II as shown below: 
n = exp | -(4.55 x ile coca a c ro al 


A table of values is not provided as this model is not applicable to 


the ocean. The values of n approach zero very rapidly. 


hs, 


Iti. DATA AND ANALYSIS 


Air Weather Service, Northern Hemisphere Historical Weather 
Maps were used to select a time interval with desirable conditions 
for this study during the period for which wave records were avail- 
able. The interval 5 January 1948 through 12 January 1948 was used 
because of a cyclonic disturbance located approximately 1500 nm 
northwest of Point Sur. This disturbance, which varied in strength, 
produced time varying wave energy. 

Microfilm reproductions of United States Weather Bureau, North 
American Surface Weather Charts were used as a source of wind-field 
data. Only actual ship wind reports were used and given a minor 
amount of smoothing. The result was a realistic wind field com- 
parable to that employed for DSA II by the French Meteorological 
service. 

Wave records were obtained from a Mark III pressure wave sensor 
installed near the bottom at Point Sur, California by the Department 
of Engineering, University of California as part of an investigation 
of surface waves conducted for the Bureau of Ships, United States 
Navy. The sensor was located in sixty-five feet of water [Isaacs 
and Saville 1949] on a bearing of 240° true from Point Sur. The 
recorder operated at a tape speed of three inches per minute for 


20 minutes out of each eight hour interval [Wiegel 1949]. 


A. HINDCAST PREPARATION 
Unlike the 'PNJ' method DSA II requires no fetch definitions 


but instead considers the action of the space- and time-varying winds 
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on each wave component. Therefore results are much less subject to 
the interpretation of the hindcaster. 
The procedure for hindcast preparation is summarized by the 


following expression: 
; dl oq] ply 
a ae it. 
oe E{¢ a 2 Pr ani ker 
0 S 


where E. is the hindcasted energy in feet square for a given period 
band at the sensor position off Point Sur; & is a summation over the 
various azimuths of energy approach to Point Sur; Gis the generation 
and decay function which, when integrated over duration, yields the 
energy in ae GOT er arriving from a specified azimuth; 2 con- 


O 
pensates for a 20 separation between azimuths instead of the 1 


specified in DSA II; C is a unit conversion (see Appendix A) from 


dl 
-L 2 
centijoules/cm to the feet squared energy representation; a7.sCts 
L 
the correction for refraction; and y— corrects for shoaling. 
S 


1. Wind-Field Plot 
Ship wind reports were transferred from microfilm to a 
plot as shown in Figure (1), which illustrates typical data density 
for the six hourly wind-field. Due to the limited detail of the 
Wind-field, azimuths from Point Sur were positioned every 20° true 
instead of every ten degrees as called for in DSA II. Azimuths were 
limited to the north by Point Reyes which bears 325° true from Point 
Sur. Bearing 320° true from Point Sur is a 300 ft shoal off Point 
Reyes which refracts waves whose period is greater than 11 seconds. 
2. Propagation Diagrams 
The propagation diagram is an extremely valuable means of 


establishing the wind history of waves arriving at an observation 


fae 


Figure 1. Sample Wind Field Plot 
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point along a particular azimuth. Figure (2) depicts a portion of 
one diagram. Actual diagrams were plotted continuously for the 
interval of 5 through 12 January 1948 to a range of 2000 nm from 
Pornt Sur. 

Each propagation diagram gives the wind conditions as a 
function of time and distance from the observation point along one 
azimuth in the wind field plot. For each increment of range (100 nm) 
and time (6 hrs) both 9 and W were plotted. Here, 9 is the angle 
between the wind direction and the direction of energy propagation 
along an azimuth and W is the wind velocity in knots. Isotachs 
were contoured every five knots of W. Wind directions were color 


coded into the following DSA II categories: 


Description 8 or W Limits 
Generating 9< Do 
Favorable 20° < @ < 60° 
Favorable-Fair 60° < @Q < 90° 
Contrary-Fair 90° < @ < 120° 

O O 
Contrary 120 <@< 180 
Calm W< 10 kts 


3. Energy Computation 


Wave record availability within each eight hour interval 
dictated hindcast times. DSA II is designed to consider waves whose 
central periods are 20, 14, 10, and 7 seconds respectively. The 7 
second waves were not used as they were obscured by hydrodynamic 
filtering and noise on the wave records. 

A propagation line overlay was prepared on the basis of 


the group velocity of the three wave bands considered (see Figure 3). 
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Figure 2. Sample Propagation Diagram. 
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Figure 3. Propagation Line Overlay 
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The advantages of the propagation diagram in combination with the 
overlay are evident. The corresponding ranges and times of each 
wave~energy packet are shown, and a hindcast can be computed for 
any given time on the diagram. 

Generation was presumed to occur when @ < 60° and one of the 


following minimum W conditions was met or exceeded: 


T (sec) Generating 0 Favorable Q 
20 40/45 kts 45 kts 
14 30/8 5% kts 35 ketes 
10 25 kts 25/30 kts 


When minimum generating or favorable winds were present 
the DSA II generation diagram for the appropriate period was entered 
with @, W, and t. Subsequent intervals with 6 < 60° either increased 
the energy density or left it unchanged depending upon W. 

When 9 was greater than 60° or W < 10 knots decay was pre- 
sumed to occur. One of the four attenuation models below was then 
employed for the appropriate conditions: 

Unmodified DSA MM fequation®(saectres 4am 

DSA II equation with 1/10th of actual duration. 

DSA II equation with 1/100th of actual duration. 

Garswell@decay function (see E12 B. 3.). 

This stepwise process was continued until the energy arrived 
at the observation point. Due to the large nearly stationary high 
pressure system to the southwest all hindcasted wave energy came from 
azimuths Dame 290°, and 310°. The reported winds were not strong 
enough at any of the hindcast times to generate waves with a period 


of 20 seconds according to the DSA II generation model. 
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The refraction coefficients (d1)/d1) below were applied to 
the energy arriving from the various azimuths [Wiegel 1964, chapter 
7]. The coefficient (d1)/d1) is the ratio of orthogonal separation 


in deep water to that at the sensor location. 


Wave Period (sec) 
Azimuth 


270° 





Within each period band the energy from the various azimuths 
was summed. The DSA II method considered azimuths of wave approach 
for every ten degrees. However in this problem the azimuths were 20" 
apart and therefore the energy in each band was doubled to obtain the 
correct total. The resulting energy density in each band was converted 
to the feet squared representation (see Appendix A). 

The arriving energy was corrected for shoaling through 
multiplication by the appropriate shoaling factor (L,/L,) below, 
where Ly is deep water wave length and L. is wave length in 65 feet 


of water (at the sensor position). The result is E. for comparison 


to the observed value of E at the sensor position. 


Period (sec) 10 14 20 


Shoaling Factor 1.38 1.68 Varo: 
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B. ANALYSIS OF WAVE RECORDS 

A grid overlay was used to sample the analogue wave records at 
two-second intervals. This provided a time series of 512 data points 
per 20 minute wave record. The data points were read in arbitrary 
units and converted to feet of head. The mean reading was computed 
and subtracted from each value, the result being a series of zero 
mean. 

Fourier amplitudes (coefficients) were obtained from the raw data 
using subroutine RHARM on the Naval Postgraduate School IBM 360 com- 
puter. RHARM is based upon the fast Fourier transform developed by 
Cooley and Tukey [1965]. The amplitudes were squared and summed into 


the following period bands: 


DSA II Wave Period Band Limits Bandwidth 
Description sec USEC ) (sec) 
Long 20 16.5-23.5 u 
Medium 14 lio ler 2 5 
Short 10 se ae ool 3 


The energy in each band is multiplied by cosH? (2rd/L. ) 
[Wiegel 1964, Appendix 1] to obtain the energy at the surface, 
where L. is the wave length at the sensor position and d is the 
depth of the pressure sensor (65 ft). The resulting energy is the 


observed energy for comparison to the hindcasted energy. 
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IV. GENERATION 


Two assumptions were made that affect the amount of energy 
generated. First, the reported wind field was presumed to represent 
the actual winds present at a constant reporting height above the 
water suitable for use with DSA II. Second, the DSA II generation 
model was presumed to be correct and was used throughout so that decay 
parameters could be varied. 

The first assumption was examined to define the effect of errors 
in the wind field. The following table indicates the approximate 
percentage error in wave energy which would result from a reported 
"senerating'" wind being 10% higher than the actual wind for a 


aquration of 10 hours. 


Actual Reported her LOU 
Wind Wind 


2a aU LNneo 


If the reported wind is considered as a mean for the six hour 
period, then speed variations during the interval would tend to make 
the hindcasted energy too small. The effect of this phenomena seems 
to be minor. 

ft is apparent that a prerequisite for this type of study is an 


accurately defined wind field. 
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V. SUMMARY OF RESULTS 


Figure (4) presents the results of the analyses done using the 
unmodified DSA II equation. The line marked "equal values" represents 
the position that should be occupied by all points if generation, decay 
and wind-field were correct and no scatter existed. This line is shown 
on each figure as a reference since various scales are used for the 
hindcast energy depending on its magnitude. 

The DSA II equation led to hindcasts that, with the exception of 
two points, were zero. It was concluded that this model provided too 
rapid an energy decay if the generation function is appropriate. 

Figure (5) shows the effect of using the DSA II decay equation 
with duration reduced to 1/10th of the actual duration. This was the 
best set of hindcasts made. Four of the five points farthest from the 
equal value line came from two hindcast times only six hours apart. 
This makes it apfteans that the reported wind-field for this area may 
have been stronger than the actual wind-field. 

Figure (6) is the final hindcast made with a form of the DSA II 
equation. One hundredth of the duration was used in the DSA II 
equation and the hindcasts are generally too large. 

A decay model slightly more rapid than the DSA II equation with 
1/10th duration would give the best results as far as general magni- 
tude is concerned. The DSA II decay showed a good balance as far as 


period is concerned. No one band of periods forms a predominant 


grouping away from the line of equal values. 
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Using DSA II Decay Equation 
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Observed Energy vs Hindcast Energy 
Using DSA II Decay Equation with 1/10th 
Actual Duration 


Y 
2 0s 
oe 
on O 
ae or 
ei, ee he, 
2 Oe ets 
=. wn 
o < 
oO 
Oo 
q 
q 
q 
oO 
. q 
o 
1.8 





Observed Energy (ftsq) 


32 


asd 


Com 


| Wea 5 


12s 


Hindcast Energy (ftsq) 


Figure 6. Observed Energy vs Hindcast Energy 
Using DSA II Decay Equation with 1/100th 


Actual Duration 
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Hindcast Energy (ftsq) 


Figure (7) presents hindcasts made with Carswell decay. Like the 
DSA II equation, decay is modeled as a function of exp [-(77)]. The 
difference is the removal of the wind dependent term, which did not 
significantly degrade the quality of the hindcast. Snodgrass et al. 
[1966, p. 493] also observed in the Pacific Ocean that swell attenu- 
ation seems to be independent of the winds encountered during 
propagation outside the generating area. 

The Groen decay function and the Phillips' scattering model were 
not used since visual inspection showed that both would make the 
hindcasts too large due to their very slow rate of decay. 

The model based on Snodgrass et al. [1966] (see II. B. 2.) was 
not used for hindcasting. In this model, n varies as exp [-(T)] and 
thus attenuates small period waves more rapidly than is observed. 

No energy was hindcast in the 20 second period band but a small 
amount was observed. It is probable that this energy arrived from 
distant regions in a manner similar to that noted by Snodgrass et al. 


[1966]. 
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VI. CONCLUSIONS 


The attenuation model which best fits the data obtained is 
n = exp (-20 r t), where n is a ratio of energy at the start of 
decay to that remaining afterwards, T is wave period in seconds, and 
t is decay duration in hours. 

The advantage of a wind dependent term, as used in DSA II, in the 
decay model is unclear from the data obtained. 

It was not possible to relate attenuation models to theories of 
decay in the ocean in this thesis. 

Future work should employ a more detailed and more precise wind 
field. Also it is possible to reduce the influence of the generation 
model through analysis. Recent innovations in computer prepared wind 
field products may be a partial answer to the wind-field problem. 
Density of reporting ships has also increased since 1948 which is an 
important gain. Knowledge of the rate of ocean wave attenuation may 


provide the key to understanding the mechanisms of decay. 
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APPENDIX A 


EQUATION OF ENERGY REPRESENTATIONS 
The French Meteorological Service representation of energy as R(X) 
1) Sanu powuseion can be converted to the amplitude squared notation 
of E(X) in fe using the following: R(X) = k = og E(X) where pg 
(sea water) is 64 ieee and k is a constant unit conversion. Intro- 


ducing appropriate unit conversions provides the following: 


x 2 
(E(x) ft? ] = on oa 2B) | R(X) een Be joule | 


centijoules 


is 0.0929 meter S 0.7376 £t-i 
a ft? joule 


E(X) £t- = 0.214 R(X) aan betes / an. 
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APPENDIX B 


TABULAR DISPLAY OF SPECTRAL COMPONENTS 
All dates refer to January 1948. Values given are energy in units 
of feet squared within the period band indicated. 
1. Observation. 
2. Hindcast using unmodified DSA II equation. 
3. Hindcast using! DSA IT equation with l710th!of actual duration, 
4. Hindcast using DSA II equation with 1/100th of actual duratilounems 


5. Hindcast using Carswell decay function. 
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® Observation. 


Period 
Date/Time 
07/1610 0.05 0.83 1.70 
07/2215 0.59 
08/0415 .0> 
08/1015 .06 


08/1615 0? 


09/1000 mus 


10/1645 02 
11/1615 .03 


12/1045 .00 
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i 


Hindcast using unmodified DSA II equation. 


Date/Time 

Long 
07/1610 O20 
07/2215 0.0 
08/0415 0.0 
08/1015 0.0 
08/1615 O20 
09/1000 0.0 
10/1645 0.0 
11/1615 0.0 
12/1045 0.0 
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Period 


Medium 


0 


0 


50) 


24 


. Ge 


3. Hindcast using DSA II equation with 1/10th of 
actual duration. 


Period 

Long Medium Short 
07/1610 0.0 0.0 0.0 
Om/ 225 0.0 1.80 0.62 
08/0415 0.0 5.49 Du? 
08/1015 0.0 13.70 5.60 
08/1615 0.0 7s 3.46 
09/1000 FO 222 4.36 
10/1645 0.0 16.17 03 
11/1615 0.0 2.49 0.61 
12/1045 0.0 2.36 0.44 
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4, Hindcast using DSA II equation with 1/100th 
of actual duration. 


Period 

Date/Time 

07/1610 O50 0.0 6.49 
07/2215 0.0 6.54 14.63 
08/0415 0.0 10.08 14,08 
08/1015 0.0 ad = 55) ero 3 
08/1615 0.0 29.05 27.90 
09/1000 0.0 142 Dorr 
10/1645 0.0 77.00 ‘aps eg 
11/1615 0.0 19.50 26.91 
12/1045 0.0 9.08 17.68 
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oY 


Hindcast using Carswell decay function. 


Date/Time — 
07/1610 0.0 
OW/ZZTD: 0.0 
08/0415 0.0 
08/1015 0.0 
08/1615 0.0 
09/1000 0.0 
10/1645 0.0 
11/1615 0.0 
12/1045 | 0.0 
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Period 


Medium 
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7.26 
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oes 


14.51 


4.51 


41.92 
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APPENDIX C 


ATTENUATION MODEL TABLES 
Carswell Attenuation Function. 
Groen Attenuation Function. 
New Zealand - Alaska Attenuation Observations. 


Phillips' Scattering Decay. 
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1. Carswell attenuation function. 
Values of n are presented as functions of period (T) and decay 
duration (t). Carswell attenuation was applied when neither favorable 


nor generating winds of greater than 10 knots were reported. 


Period 
t (hours) 
Wy 2 G99 04. 96 0.96 
L GE76 0.96 Os 
2 0'.. 96 0.95 0.86 
5 O.9F OvmZ 0.68 
10 0.82 0.67 0.47 
20 0.68 0.46 O22 
30 0.5/7 0.38 0.10 
50 0.89 0.14 0.02 
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2. Groen attenuation function. 
Values of n are presented as functions of period (T) and decay 
duration (t). This function was not used since it was apparent by 


inspection that the resulting hindcasts would be much too large. 


Period 

Ea(houre® Long Medium Short 
1/2 1.00 1.00 1 (O10 

1 iL Ae) 0 1.00 1.00 

9 1.00 1.00 1.00 

5 1.00 1.00 1.00 

10 1.00 0.99 0.99 

20 0.99 0.99 0.98 

30 0.99 0.98 0.97 


50 O98 Omg, 0296 
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3. New Zealand - Alaska attenuation observations. 

Values of n are presented as functions of period (T) and decay 
duration (t). This function is slightly more rapid in its attenuation 
of energy than the Carswell model. Compared to the results of this 
study decay is of the right order of magnitude. However, the frequency 
distribution does not appear correct. This study shows that short waves 


attenuate more rapidly than medium waves rather than as shown below. 


Period 
t (hours) 
1/2 0.9259 0.9465 
I! 0.8573 0.8958 
2 0.7349 0.8025 
D 0.4630 0.5770 
10 0.2144 0,.3329 
20 0.0460 0.1108 
30 0 10099 0.0369 
50 0.0004 0.0041 
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4. Phillips Scattering Decay 
Values of n are presented as functions of period (T) and decay 
duration (t). This function was not used for hindcasts since, by 


inspection, the rate of attenuation it models is much slower than that 


actually observed. 
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